Calculation of nuclear-spin-dependent parity nonconservation in s-d transitions of 

Ba + , Yb + and Ra + ions. 
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We use correlation potential and many-body perturbation theory techniques to calculate spin- 
independent and nuclear spin-dependent parts of the parity nonconserving amplitudes of the tran- 
sitions between the 6S1/2 ground state and the 5^3/2 excited state of Ba + and Yb + and between 
the 7si/2 ground state and the 6^3/2 excited state of Ra + . The results are presented in a form con- 
venient for extracting of the constants of nuclear-spin-dependent interaction (such as, e.g., anapole 
moment) from the measurements. 
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I. INTRODUCTION 

The study of the parity nonconservation (PNC) in 
atoms is a low-energy, relatively inexpensive alternative 
to high-energy search for new physics beyond the stan- 
dard model (see, e.g. [l|). The most significant recent 
achievement on this path is the very precise measure- 
ments of the PNC in cesium Q ■ The cesium PNC exper- 
iment together with its interpretation in terms of 
nuclear weak charge provides the best current atomic test 
of the standard model (see also review @). It is also the 
only measurement of the nuclear anapole moment which 
is produced by the PNC nuclear forces Q. The extrac- 
tion of the weak nuclear charge from the PNC measure- 
ments relies on atomic calculations. Cesium atom has 
the simplest electron structure among all heavy atoms 
which were used or considered for the PNC measure- 
ments. Still it took considerable efforts of several groups 
of theorists to bring the accuracy of the calculations in 
line with the accuracy of measurements and provide reli- 
able interpretation of the measurements in terms of the 
standard model and possible new physics beyond it 0- 
5]. It is widely believed now that it would be hard to 
compete with cesium experiment in terms of accuracy 
of interpretation of the PNC measurements. Therefore, 
the study of PNC in atoms is mostly focused now in two 
directions: (i) the measurements of the PNC ratio for 
a chain of isotopes which was first proposed in Ref. Q, 
and (ii) the measurements of the nuclear-spin-dependent 
PNC, like e.g. the contribution from nuclear anapole mo- 
ment (see, e.g. reviews .6, 9]). The study of the PNC for 
a chain of isotopes does not require atomic calculations 
and can deliver useful information about either neutron 
distribution or new physics beyond standard model (see, 
e.g. (Tol - [l^ |). The measurements of anapole moment does 
require atomic calculations but high accuracy is not crit- 
ical here. 

Ba + , Yb + and Ra + ions considered in present paper 
are good candidates for both types of the experimen- 
tal studies. Ba and Yb both have seven stable isotopes 
with large difference in neutron numbers AN max = 8. 



Radium has several long-living isotopes. There are two 
stable isotopes for each of the Ba and Yb atoms ( 135 Ba, 
139 Ba, 171 Yb and 173 Yb) which have non-zero nuclear 
spin. There are also isotopes of Ra with non-zero nuclear 
spin ( 223 Ra, 225 Ra, 229 Ra). In all cases nuclear spin is 
provided by valence neutron. This is especially interest- 
ing since it allows one to measure the strength of the 
neutron-nucleus PNC potential Q (the anapole moment 
has been measured only for the 133 Cs nucleus which has 
valence proton). 

Finally, Ba + and Ra + ions have electron structure sim- 
ilar to those of cesium atom. This means that the accu- 
racy of the interpretation of the PNC measurements can 
be on the same level as for cesium. Moreover, it can 
be further improved with the use of the experimental 
data [3. 

The use of Ba + in the PNC measurements was first 
suggested by Fortso n 1141 . The work is in progress at 
Seattle (see, e.g. fl5lll6l]) but no PNC results have been 
reported yet. Similar approach is now considered for 
the measurements of PNC in Ra + ion at KVI [l8|, . 
It is important that in Ra + the PNC effects are about 
20 times larger than in Ba + . There are plans to mea- 
sure PNC in Yb + at Los Alamos [loj]. Note that the 
PNC measurements for neutral ytterbium are in progress 
at Berkeley and first PNC results were recently re- 
ported [ll|. The PNC measurements for the Yb + ion 
would provide an important consistency test for the mea- 
surements and their interpretation. 

Calculations of the spin-independent PNC amplitude 
for Ba + and Ra + were performed in our early work (l3l | 
and in [17|. Calculations for Ra + were later performed in 
(18TJ and [22[ . The only calculation of the spin-dependent 
PNC in Ra + was recently reported by Sahoo et al 
To the best of our knowledge, no PNC calculations for 
Yb + have been published so far. 

In present paper we calculate both spin-independent 
and spin-dependent PNC amplitudes simultaneously us- 
ing the same procedure and the same wave functions. In 
this approach the relative sign of the amplitudes is fixed. 
This allows for unambiguous determination of the sign 
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of the spin-dependent contribution. The constant of the 
spin-dependent interaction can be expressed via the ratio 
of the two amplitudes. This brings an extra advantage of 
more accurate interpretation of the measurements. The 
accuracy of the calculations for the ratio of the PNC am- 
plitudes is usually higher than that for each of the am- 
plitudes. This is because the amplitudes are often very 
similar in structure and most of the theoretical uncer- 
tainty cancels out in the ratio. 

Since we focus on the calculation of the nuclear-spin- 
dependent PNC amplitudes where high accuracy of cal- 
culations is not needed, we don't include some small cor- 
rections, like some classes of diagrams for higher-order 
correlations, Breit and quantum electrodynamic (QED) 
corrections, etc. Instead, we make sure that all leading 
contributions are included exactly the same way for both 
spin- independent and spin-dependent PNC amplitudes 
which is important for the cancelation of the uncertainty 
in the ratio. 



II. THEORY 

Hamiltonian describing parity-nonconserving electron- 
nuclear interaction can be written as a sum of spin- 
independent (SI) and spin-dependent (SD) parts (we use 
atomic units: ft = |e| = m e = 1): 
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where Gf ~ 2.2225 x 10 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 
er 







and 75 are the Dirac matrices, I is the nu- 



clear spin, and p(r) is the nuclear density normalized to 
1. The strength of the spin-dependent PNC interaction 
is proportional to the dimensionless constant k which 
is to be found from the measurements. There are three 
major contributions to h arising from (i) electromagnetic 
interaction of atomic electrons with nuclear anapole mo- 
ment 24], (ii) electron- nucleus spin-dependent weak in- 
teraction, and (iii) combined effect of spin-independent 
weak interaction and magnetic hyperfine interaction [25[ 
(see, also review [6j). In this work we do not distinguish 
between different contributions to h and present the re- 
sults in terms of total >c which is the sum of all possible 
contributions. 

Within the standard model the weak nuclear charge 
Qw is given by [26j 
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-0.98777V + 0.0716Z. 



(2) 



Here N is the number of neutrons, Z is the number of 
protons. 

The PNC amplitude of an electric dipole transition 



between states of the same parity \i) and |/) is equal to: 
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where d—~e J^i r i is the electric dipole operator, \a) = 
\J a F a M a ) and F = I+J is the total angular momentum. 

Applying the Wigner-Eckart theorem we can express 
the amplitudes via reduced matrix elements 
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Detailed expressions for the reduced matrix elements of 
the SI and SD PNC amplitudes can be found e.g. in 
Refs. [27] and [28]. For the SI amplitude we have 
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For the SD PNC amplitude we have 
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For the case of the 5c? — 6s transitions considered in 
present paper (or 6d — 7s in the case of Ra + ) it is conve- 
nient to break expression ((5]) into four parts: 

(5d 3/2 ,F f \\dsD\\6a,F i ) = S 1 + S2 + S 3 + S i , (7) 
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where 
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Here c m (Ff,Fi) (m = 1,2,3,4) are coefficients which can 
be reconstructed using ([5]). The terms Si, S 2 , S 3 , S4 differ 
by the order of the operators d and ap and by the states 
in the summation which are either npi/ 2 or np 3 / 2 states. 
To know the relative values of these terms is important 
for the analysis of the accuracy of the calculations. 



TABLE I: Ionization energies of lowest s,p and d states of 
Ba + , Yb + and Ra + in different approximations (cm -1 ). 



Ion 


State 


RHF 


Brueckner 


Experiment [36] 


Ba + 


6«l/2 


75340 


80815 


80687 




6P1/2 


57266 


60571 


60425 




6P3/2 


55873 


58848 


58735 




5d 3/2 


68139 


76318 


75813 


1 D 


DSl/2 


yU /oy 


GQ/1 11 

yy4/ f 


yozU f 




6P1/2 


66087 


70728 


71145 




6P3/2 


63276 


67101 


67815 




5d 3/2 


66517 


75551 


75246 


Ra+ 


7S1/2 


75898 


82032 


81842 




7P1/2 


56878 


60715 


60491 




7P3/2 


52906 


55753 


55633 




6^3/2 


62356 


70091 


69758 


TABLE II: Rescaling factors for the correlation potential E. 


Ion 


Sl/2 




P1/2 


P-i/2 d 3 / 2 


Ba+ 


0.978 




0.960 


0.964 0.941 


Yb+ 


0.862 




1.081 


1.170 0.968 


Ra+ 


0.970 




0.946 


0.960 0.959 



III. CALCULATIONS 

To perform the calculations we follow an ah initio ap- 
proach which uses the correlation potential method [29] 
and the technique to include higher-order correlations de- 
veloped in Refs. [30l - l32l |. 

Calculations start from the relativistic Hartree-Fock 
(RHF) method in the V N ~ Y approximation. This means 
that the initial RHF procedure is done for a closed-shell 
atomic core with the valence electron removed. After 
that, the states of the external electron are calculated 
in the field of the frozen core. Correlations are included 
by means of the correlation potential method (29|. For 
Ba + and Ra + we use the all-order correlation poten- 
tial which includes two classes of the higher-order 
terms: screening of the Coulomb interaction and hole- 
particle interaction (see, e.g. [3(| for details). For Yb + 
we use the second-order correlation potential E^ 2 ). The 
reason for different approaches is due to different elec- 
tron structures of the ions. The all-order technique de- 
veloped in [30Tj32l ] works very well for alkali atoms and 
similar ions in which the valence electron is far from the 
atomic core and higher-order correlations are dominated 
by screening of the core- valence residual Coulomb inter- 
action by the core electrons. For atoms and ions simi- 
lar to Yb + , in which an external electron is close to the 
core and strongly interacts with its electrons, a different 
higher-order effect described by the ladder diagrams [33| 
becomes important. The applicability of the technique of 
Rcf. [33] to Yb + needs further investigation. Meanwhile, 
the use of the second-order X^ 2 ) leads to sufficiently good 
results. Note that an external electron in Ba + and Ra + 



ions is also closer to atomic core than in neutral alkali 
atoms Cs and Fr. This means that inclusion of ladder 
diagrams might be a way to improve the accuracy of cal- 
culations for the ions as well. This question also needs 
further investigation. 

To calculate £ or £( 2 ))we need a complete set 

of the single-electron orbitals. We use the B-spline tech- 
nique [3 l] to construct the basis. The orbitals are built as 
linear combinations of 50 B-splines of order 9 in a cavity 
of radius 40as. The coefficients are chosen from the con- 
dition that the orbitals are the eigenstates of the RHF 
Hamiltonian Hq of the closed-shell core. The second- 
order operator X^ 2 ) is calculated via direct summation 
over B-spline basis states. The all-order £(°°) is calcu- 
lated with the technique which combines solving equa- 
tions for the Green functions (for the direct diagram) 
with the summation over complete set of states (exchange 
diagram) [30l ] . 

The correlation potential X is then used to build a 
new set of single-electron states, the so-called Brueck- 
ner orbitals. This set is to be used in the summation in 
equations (0, and ((6|). Here again we use the B-spline 
technique to build the basis. The procedure is very sim- 
ilar to constructing of the RHF B-spline basis. The only 
difference is that new orbitals are now the eigenstates of 
the Hq + S Hamiltonian. 

Brueckner orbitals which correspond to the lowest va- 
lence states are good approximations to the real physical 
states. Their quality can be tested by comparing ex- 
perimental and theoretical energies. The energies of the 
lowest states of Ba + , Yb + and Ra + in RHF and Brueck- 
ner approximations are presented in Table |U One can 
see that inclusion of the correlations leads to significant 
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improvement of the accuracy in all cases. The devia- 
tion of the theory from experiment is just fraction of a 
per cent in the case of Ba + and Ra + where an all-order 
i s U sed and does not exceed 1.3% for Yb + where 
the second-order E^ 2 ) is used. 

The quality of the Brueckner orbitals can be further 
improved by rescaling the correlation potential E to fit 
the experimental energies exactly. We do this by replac- 
ing the Ho + E with the Ho + AE Hamiltonian in which 
the rescaling parameter A is chosen for each partial wave 
to fit the energy of the first valence state. The values 
of A are presented in Table [TTJ Note that these values 
are very close to unity. This means that even without 
rescaling the accuracy is good and only a small adjust- 
ment of the value of E is needed. Note also that since 
the rescaling procedure affects not only energies but also 
the wave functions, it usually leads to improved values 
of the matrix elements of external fields. In fact, this is 
a semi-empirical method to include omitted higher-order 
correlation corrections. 

Matrix elements of the Hsi, Hgu and electric dipole 
operators are found by means of the time-dependent 
Hartree-Fock (TDHF) method [H, [H| extended to 
Brueckner orbitals. This method incorporates to the 
well-known random-phase approximation (RPA) dia- 
grams including exchange. In the TDHF method, the 
single-electron wave functions are presented in the form 
ip = tpo + Sip, where 4*0 is the unperturbed wave func- 
tion. It is an cigcnstate of the RHF Hamiltonian Ho: 
(Ho — £o)V'o = 0. Sip is the correction due to external 
field. It can be found be solving the TDHF equation 

{H - e )5*l> = -6eip ~ Fip - SV^ 1 ^, (12) 

where Se is the correction to the energy due to external 
field (Se = for all above mentioned operators but it 
is not zero for the hyperfine interaction which we will 
need for the analysis of accuracy), F is the operator of 
the external field, and SV 1 * -1 is the correction to the 
self-consistent potential of the core due to external field. 

The TDHF equations are solved self-consistently for 
all states in the core. Then the matrix elements between 
any (core or valence) states n and m are given by 

(i> n \F + 5V N - l \^ m ). (13) 

The best results are achieved when ip n and ip m are 
the Brueckner orbitals computed with rescaled correla- 
tion potential E. 

We use equation (fT3|) for all weak and electric dipole 
matrix elements in evaluating the SI and SD PNC am- 
plitudes © and ©. 

IV. ACCURACY OF CALCULATIONS 

The accuracy of the results obtained via direct sum- 
mation over physical states with the use of expressions 
like ([3]) is determined by the accuracy for the energies, 



TABLE III: Electric dipole matrix elements. Comparison of 
present calculations with experiment or most complete other 
calculations. 





TV ariciTi An 
LI (AIIOILILJII 


L Ilia WUIri 


Other 


Ba+ 


6s 1/2 - 6pi /2 
6s 1/2 - 6p 3 / 2 

5^3/2 - 6pi/ 2 
5^3/2 - 6p 3 /2 


3.32 
4.69 
3.06 
1.34 


3.36(4)" 
4.55(10)" 
2,.l<±($) b 
1.54(19)" 


Yb+ 


6s 1/2 - 6pi /2 
6s 1/2 - 6p 3/2 
5d 3/2 - 6p 1/2 
5d 3/2 - 6p 3/2 


2.72 
3.84 
3.09 
1.36 


2.471(3) c 
3.36(2) d 
2.97(4) c 
1.31 e 


Ra+ 


7Sl/2 - TPl/2 
7Si/2 - 7p 3 /2 

6d 3 / 2 - 7pi/ 2 
6d 3 /2 - 7p 3/2 


3.24 
4.49 
3.56 
1.51 


3.254' 
4.511 / 
3.566 / 
1.512' 


"Experiment, Ref. [37|. b Experiment, Ref. [la]. 
Experiment, Ref. [3lj3i]™Experiment i Ref. [S|. 'Theory, 
RefTg^^heory, Ref. g^. 


TABLE IV: 
Comparison 


Magnetic dipole hyperfine constants A (MHz), 
of present calculations with experiment. 


Ion 


State 


This work 


Experiment 


13b Ba+ 


6Sl/2 
6pi/2 
6P3/2 
5^3/2 


3671 
668 
131 
161 


3593.3(2.2)" 
664.6(0.3) 6 
113.0(0.1)' 
169.5892(9) c 


m Yb + 


6S1/2 
6P1/2 
6P3/2 
5^3/2 


13217 
2533 
388 
291 


12645(2)" 
2104.9(1.3) d 
877(20) c 
430(43/ 


^ 3 Ra+ 


7S1/2 
7P1/2 
7P3/2 
6^3/2 


3537 
679 
69.8 
57.8 


3404(2) 9 
667(2) 9 
56.5(8) s 
77.6(8) ft 



"Ref. EJ b RefJipl <Ref. 0. d Ref. [Hf "Ref. gfj. 
'Ref. H^ 9 Ref. 5aBl- ft Rescaled from 213 Ra [TJ] using 
magnetic moments from [50l |. 



electric dipole and weak matrix elements. We start from 
the notion that for the PNC amplitudes considered in 
present work the summation over intermediate p-states 
is strongly dominated by the Qpi/2 and 6p 3 / 2 states for 
Ba + and Yb + and by lp\/2 and 7p 3 / 2 states for Ra + . 
Corresponding contributions constitute 70 to 90% of the 
total PNC amplitude. Therefore, it is sufficient to com- 
pare with experiment energies and matrix elements in- 
volving these p-states. The energies and electric dipole 
matrix elements can be directly compared with experi- 
ment while standard practice of comparing experimental 
and theoretical hyperfine structure can be used to test 
the accuracy of the weak matrix elements. 

To improve the accuracy for the amplitudes the ener- 
gies of the 6s, 6pi/2, 6P3/2 and 5(^/2 states (7s, 7p\/2, 
7p3/2 and 6^3/2 for Ra + ) are fitted exactly in our calcu- 
lations using rescaling of the correlation potential E as it 
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has been described in previous section. 

Calculated and experimental El-transition amplitudes 
are presented in Table Hill Note that we need comparison 
with experiment only for estimation of the accuracy of 
our calculations. Therefore, a comprehensive review of 
the experimental and theoretical data available for the 
ions goes beyond the scope of present work. We only 
compare our results with the most accurate experimental 
data or with the most complete other calculations where 
the experimental data are not available. Good reviews of 
the electric dipole transition data in Ba + and Yb + can 
be found in Ref. PJ] and HI. 

The data in Table Mil shows good agreement between 
theory and experiment for most of the amplitudes, al- 
though the accuracy for the amplitudes involving the Pz/2 
states is lower than that for the P1/2 states. 

Table IIVI shows theoretical and experimental data on 
the hyperfme structure constants of the low states of 
Ba + , Yb + and Ra + . Here again we only compare our cal- 
culations with the most accurate experimental data. A 
review of the available experimental and theoretical data 
for Ba + can be found in Ref. 51] . The data in Table IIVI 



shows several trends: (i) the accuracy is good for S1/2 and 
P1/2 states, especially in the cases of Ba + and Ra + , (ii) 
the accuracy for Yb + is lower than that for Ba + and Ra + , 
(iii) the accuracy for p 3 / 2 and d 3 / 2 states is lower than 
that for the S1/2 and P1/2 states. The largest discrepancy 
is for the hfs of the 63/2 state of Yb + where theory and 
experiment differ almost three times. Note that the most 
complete calculations of Ref. [4l| give the result which is 
close to our theoretical value rather than to the experi- 
ment. In principle, the discrepancy can be explained by 
configuration mixing involving configurations with exci- 
tations from the 4/ subshell. Neither our present calcula- 
tions nor those of Ref. (4l| include this mixing explicitly. 
The configuration interaction calculations based on tech- 
nique developed in Ref. [H, [HI which treats Yb + as a 
system with fifteen valence electrons show that the hfs of 
the 6P3/2 state is indeed very sensitive to the configura- 
tion mixing. One can find such mixing which reproduces 
the experimental hfs exactly while the accuracy for the 
energy and for the g- factor of the 6p 3 / 2 state is also good. 
However, the results are inconclusive due to strong insta- 
bility of the hfs of the 6P3/2 state. We can only say that 
the configuration mixing can explain current experimen- 
tal value of the hfs of 6p 3 / 2 state but we cannot prove 
that this explanation is correct. Since the disagreement 
between theory and experiment for the hfs of the 6^3/2 
state of Yb + is the main factor contributing to the un- 
certainty of the calculations for Yb + , it would be useful 
to remeasure the hfs of this state. 

The fact that the accuracy for the P1/2 and p 3 / 2 states 
is different complicates the analysis of the accuracy for 
the PNC amplitudes. There is cancelation between terms 
containing matrix elements with the P1/2 and pz/2 states. 
In the end of section |TT] we introduced the notations 
51,5*2,53 and 5 4 for these terms (see Eqs. (|8l9ll0llljl ). 
The terms involving the P1/2 states are 5i and 5 3 , the 



terms with the p^/2 states are 52 and S4. Table IVl shows 
the Si, 52, 53, 54 contributions to the reduced matrix el- 
ements of the nuclear-spin-dependent PNC interaction in 
some hfs components of the transitions in Ba + , of Yb + 
and of Ra + . One can see that the 52 term is usually 
small while the 54 term is not small. For example, for 
Yb + the contribution of the 54 term is more than a half 
of the total sum. It is clear that the accuracy of the cal- 
culations in this case will be mostly determined by the 
accuracy of the 54 term. 

To analyse the accuracy of the PNC calculations we 
need a procedure which takes into account the deviation 
of the experimental and theoretical data for the electric 
dipole matrix elements and for the hyperfme structure as 
well as the effect of partial cancelation between different 
contributions to the PNC amplitude. We do this by com- 
paring the ab initio calculations with the calculations in 
which the electric dipole and weak matrix elements are 
rescalcd to fit the experimental data. For example, as- 
suming that the weak matrix elements between two states 
are proportional to the square root of the hfs constants 
for these states we rescale them as following 



(n -ffpNC|?7l) re scaled 



A CX P A CX P 
^th^th 



;n|JT P Nc|m). (14) 



Here A cxp and are experimental and theoretical val- 
ues of the hfs constants from Table IIVI This means that 
we perform accurate rescaling for matrix elements involv- 
ing 6p!/ 2 and 6p 3 / 2 states (7p\/2 and 7p 3 / 2 for Ra + ). As 
it was stated above, this corresponds to 70 to 90% of the 
total PNC amplitude. We use the same rescaling for all 
matrix elements involving higher p states. Electric dipole 
matrix elements are also rescaled to fit the experimental 
data for the transitions between lowest states. The differ- 
ence between PNC amplitudes obtained in the ab initio 
calculations and calculations with rescaling serves as an 
estimation of the uncertainty of the calculations. 

Note that the accuracy for the relative contribution of 
the nuclear-spin-dependent interaction can be higher that 
for each of the amplitudes (see also Ref. [54| ) . As we will 
see in the next section, this is usually the case when the 
52 and 53 contributions arc both small. This is because 
these terms are exactly zero for the spin-independent 
PNC amplitudes. Therefore, the spin-dependent PNC 
amplitudes in which the 52 and 53 terms are small, are 
similar to the spin-independent amplitudes. They both 
change under scaling at the same rate which cancels out 
in the ratio. 
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TABLE V: Contributions to the reduced matrix elements 
{5d 3/ 2,-F , i||i5 r sDPNcl|6si/2,-F2> of the spin-dependent parity- 
nonconserving s-d transitions. See text for explanation of 
notations. Units: 10 -11 >«eao. 



Ion 


Ft 


F 2 


Si 


52 


S3 


S4 


Sum 







1 


0.134 


0.002 


0.000 


-0.027 


0.108 




1 


1 


-0.211 


-0.001 


0.013 


0.032 


-0.168 




1 


2 


-0.057 


0.003 


0.029 


-0.014 


-0.040 




2 


1 


0.211 


-0.002 


-0.038 


-0.009 


0.162 




2 


2 


0.127 


-0.003 


-0.038 


0.009 


0.094 




3 


2 


-0.212 


-0.002 


0.000 


0.043 


-0.171 


lvl Yb + 


1 





0.780 


0.000 


-0.306 


-0.164 


0.310 




1 


1 


0.184 


-0.008 


-0.432 


0.116 


-0.140 




2 


1 


-0.411 


-0.004 


0.000 


0.156 


-0.259 


^ y Ra+ 


1 


2 


2.021 


0.031 


0.000 


-0.119 


1.933 




2 


2 


-2.301 


-0.005 


0.265 


0.084 


-1.957 




2 


3 


-0.878 


0.044 


0.496 


-0.045 


-0.384 




3 


2 


2.058 


-0.037 


-0.593 


-0.006 


1.423 




3 


3 


1.643 


-0.036 


-0.530 


0.006 


1.084 




4 


3 


-2.500 


-0.039 


0.000 


0.148 


-2.391 



V. RESULTS 

The results of the calculations for the spin- independent 
part of the PNC amplitudes (z-components) are 

I PNC/ 



Ba^ 



Yb+ 



Ra^ 



El^(5d 3/2 - 6s) = 
0.29(2) x lQ- 12 Q w iea Ql (15) 
£l PNC (5d 3 /2 - 6a) = 
0.62(20) x l0- 12 Q w iea„, (16) 
£l PNC (6d 3/2 -7 S ) = 
3.4(1) x 10~ 12 Q w iea . (17) 



The uncertainties are estimated by comparing ab initio 
calculations with the calculations in which matrix ele- 
ments were rescaled as it was described in previous sec- 
tion. The expressions (|15I16I17[) are valid for any iso- 
topes. All dependence on nuclear number A is via weak 
nuclear charge Qw (see, @) while dependence on nuclear 
radius is negligible. To be precise, the dependence of the 
PNC amplitudes on the nuclear radius can be included 
with the help of an additional factor 



to the spin-independent amplitudes. This has two impor- 
tant advantages [54] : (i) extraction of the value of x from 
experimental data can lead to no confusion over its sign, 
(ii) the uncertainty for the value of the ratio of the spin- 
dependent and spin-independent amplitudes R is usually 
lower than for each of the amplitudes. This is because 
the two amplitudes are very similar and numerical un- 
certainty cancels out in the ratio (see also Ref. [54|). 

The total PNC amplitudes for different hfs transitions 
in Ba + , Yb + and Ra + are presented in Table IVD The 
table includes all stable isotopes of Ba and Yb which 
have non-zero nuclear spin and the most stable isotopes 
of Ra with non-zero nuclear spin. The results for other 
isotopes can be obtained by rescaling appropriate PNC 
amplitude (with required values of F± , F2 and I) using 
corresponding weak nuclear charges: 



El PNC (A 2 ) FlF2l 



(20) 



P(A 



i)F 1 F 2 r i 



!W 



(A 2 



(A, 



1 + R(Ai)f 1 f 2 i 



!W 



(^l) 



Qw(A 2 )' 



where P(Ai)p 1 p 2 j and R{A\)f 1 f 2 i are taken from Ta- 
ble [VI] and Qw(Ai) and Qw(A 2 ) are calculated using 
We stress ones more that the dependence of the 
amplitudes on the nuclear radius is much smaller than 
current theoretical uncertainty. 

Numerical uncertainties for P and R are presented in 
parentheses in Table I VII One can see that for some 
hyperfme transitions the uncertainty for R is very low. 
Comparing the data in Tables IVT1 and IVl reveals that low 
uncertainty in R corresponds to the cases when the spin- 
dependent PNC amplitude is strongly dominated by the 
sum 5*i + S4 while the sum of two other terms (S 2 and 
S3) is small. This is because strong domination of S\ + S4 
makes the spin-dependent PNC amplitude to be very 
similar to the spin-independent one where S 2 = and 
S3 = 0. In this case the rescaling changes both ampli- 
tudes at the same rate and the change cancels out in the 
ratio R. The hfs transitions with low uncertainty in R 
are good candidates for the measurements when the aim 
is extraction of h. 



A. Comparison with other calculations 



El 



PNC 



(A 2 ) = 



A x 



El™ c (Ai). 



(18) 



For cases considered in this work the maximum value of 
the correction is 0.4% (between 223 Ra and 229 Ra). For 
other cases the correction is even smaller. This is beyond 
the accuracy of present calculations. 

It is convenient to present the total PNC amplitude 
(including the spin-dependent part) in a form 



E1 PNC = P{1 + R> ( ), 



(19) 



where P is the spin-independent part (including weak nu- 
clear charge Qw) an d R is the ratio of the spin-dependent 



Table IVlTl summarizes present and past calculations of 
the spin-independent PNC s-d amplitudes in Ba + and 
Ra + . We present the results in a form of the coeffi- 
cients before weak nuclear charge Qw- These coefficients 
are practically isotope-independent. This is because the 
isotope-dependence of the PNC amplitudes is strongly 
dominated by weak nuclear charge while the dependence 
of the PNC amplitudes on the details of nuclear density 
is very weak and can be neglected on the present level of 
accuracy. 

The technique used in the present work is very sim- 
ilar to the sum-over-states approach of our previous 
paper [13j. As expected, the results are very close 
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TABLE VI: 

\5d 3/2y F 1 ) 



PNC amplitudes (z-components) 
^1/2, F2) transitions in 135 Ba + , 



for the 
137 Ba+, 

171 Yb+ and 173 Yb+ and |6d 3/2 ,Fi) -» \7s 1/2 ,F 2 ) transitions 



TABLE VIII: Reduced matrix of the spin-dependent parity- 
nonconserving s-d amplitudes in Ba + and Ra + . Units: 



in 223 Ra+, 225 Ra+ and 229 Ra+. Units: 10 



te«o. 



Ion 


Qw 


I 


Fi 


F, 


PNC amplitude 


i3b Ba+ 


-74.11 


1.5 





1 


—0.152(9) 


X 


[" -1 1 aa,iaa/a\ 1 

[1 + 0.0409(2) 3t\ 








1 


1 


—0.170(11) 


X 


1" 1 1 aa,iaa/a\ 1 

[1 + 0.0400(2)>«-] 








1 


2 


—0.059(4) 


X 


[1 — 0.021(2)>rJ 








2 


1 


0.132(9) 


X 


[1 + 0.039(l);nr] 








2 


2 


—0.152(9) 


X 


[1 — 0.023(l)>rJ 








3 


2 


\J.±o2{\3) 


X 


[1 — U.U24o(l 


137 Ba+ 


-76.09 


1.5 





1 


—0.156(10) 


X 


[1 + 0.0398(2)x-j 








1 


1 


—0.175(11) 


X 


[1 + 0.0392(3)x-] 








1 


2 


—0.061(4) 


X 


[1 — \J.U2i{2)x\ 








2 


1 


0.135(8) 


X 


r-i i a aoo / 1 \ i 

[1 + 0.038(l)>rJ 








2 


2 


a 1 tr / 1 A \ 

—0.156(10) 


X 


[1 — \J.U22{i)x\ 








3 


2 


0.1OO(10J 


X 


M n noon/ 1 ^ ,. ,1 
[1 — \J.\j2d\){± )k\ 


171 Yb+ 


-94.86 


0.5 


1 





0.59(19) 


X 


[1 + 0.030(16)>f] 








1 


1 


—0.29(9) 


X 


[1 + 0.019(2)>«rJ 








2 


1 


0.51(16) 


X 


[1 — 0.016(6)>ir] 


173 Yb+ 


-96.84 


2.5 


1 


2 


-0.41(13) 


X 


[1 + 0.022(9);nr] 








2 


2 


—0.53(17) 


X 


[1 + 0.015(8)>rJ 








2 


3 


—0.17(6) 


X 


Ti i A AAA ( C\\ 1 

[1 + 0.009(2)>«rJ 








3 


2 


0.28(9) 


X 


[1 + 0.005(3)>ir| 








3 


3 


—0.48(5) 


X 


[1 — \SX^)2{\)k\ 








4 


3 


0.3/(12] 


X 


[1 — {j.{j±0{Q)>C\ 


223 Ra+ 


-127.2 


1.5 





1 


— 3.04(9) 


X 


[1 + 0.0252(l)x-J 








1 


1 


a a / 1 a\ 

—3.40(10) 


X 


[1 + 0.0233(3)x-] 








1 


2 


— 1.18(4) 


X 


[1 — 0.0053(5)x-] 








2 


1 


2.64(8) 


X 


[1 + 0.0193(3)x] 








2 


2 


-3.04(9) 


X 


[1 - 0.0093(4)*:] 








3 


2 


3.04(9) 


X 


[1 - 0.0151(1)^] 


225 Ra+ 


-129.2 


0.5 


1 





4.37(13) 


X 


[1 + 0.0389(3)^] 








1 


1 


-2.19(6) 


X 


[1 - 0.0033(6)^] 








2 


1 


3.79(11) 


X 


[1 - 0.0149(1);*] 


229 Ra+ 


-133.1 


2.5 


1 


2 


-3.02(9) 


X 


f~l _l_ n A9A9/1 A i^l 








2 


2 


-3.97(12) 


X 


[1 + 0.0180(1)*] 








2 


3 


-1.27(4) 


X 


[1 - 0.0066(4)*] 








3 


2 


2.12(6) 


X 


[1 + 0.0146(3)*] 








3 


3 


-3.56(10) 


X 


[1 - 0.0100(3)*] 








4 


3 


2.76(8) 


X 


[1 - 0.0145(1)*] 


TABLE 


VII: 


Spi 


n-independent part 


of the parity- 


nonconserving 


s-d 


amplitudes in Ba + , 


Yb+ and Ra+. 


Units: 10 12 Qwieao- 










Ion 


Transition 




This work 




Other 


Ba+ 


5^3/2 


-6s 


1/2 




0.29(2) 


0.29", 0.304" 


Yb+ 


5rf 3 / 2 


- 6s 


1/2 




0.62(20) 






Ra+ 


6^3/2 


- 7s 


1/2 




3.4(1) 


3.3 a , 3.36 c , 3.33 d 




a Ref. [Q. ' 


Ref. 


[17(. c Ref. [l8j. 


d Ref. 22]. 



too. There is also good agreement with Sahoo et al for 
Ba + [13] and with Wansbeek et al for Ra + PH and with 
recent calculations by Pal et al [13] for Ra + . 

Tabic [VlTTI compares our calculated reduced matrix el- 
ements of the spin-dependent PNC amplitudes with the 
results of the recent calculations by Sahoo et al [23| . To 



10 12 xieao- 



Ion 


Transition 


I 


Fi 


F 2 


This work 


Ref. [23] 


13b Ba+ 


6S1/2 — 5(^3/2 


1.5 


2 


3 


-1.71 


-1.94 






1 


2 


1.62 


1.79 


139 Ba+ 


6si/ 2 - 5d 3 / 2 


3.5 


3 


3 


-1.86 


-2.07 






3 


2 


1.86 


2.11 


22b Ra+ 


7si/ 2 - 6d 3 / 2 


0.5 


1 


2 


-17.8 


-19.8 


223 Ra+ 


7s 1/2 - 6d 3 / 2 


1.5 


2 


3 


-21.1 


-23.5 






1 


2 


16.1 


20.3 


229 Ra+ 


7si/ 2 - 6d 3/2 


2.5 


2 


3 


-3.8 


-6.5 






2 


2 


-19.6 


-22.9 



make the comparison easy we have multiplied all matrix 
elements from [23] by 2 and have changed their signs. 
The former is to take into account different definition 
of *, the latter is due to the fact that we also have an 
opposite sign for the spin-independent PNC amplitude 

[Hj]. The to- 



compared to what is presented in [lSj and 
tal sign of an amplitude is not fixed and can be changed 
arbitrarily. Note however that the relative sign of the SI 
and SD PNC amplitudes is not arbitrary and the sign 
can only be changed for both parts of the amplitudes 
simultaneously. 

Given that the accuracy of the present calculations is 
few per cents and similar accuracy should be expected for 
[23j the results presented in Table [Villi are in a reasonable 
agreement with each other. Comparison of the data in 
Table IVIlTl and TablelVlshows that the difference between 
our results and those of Sahoo et al is larger for cases 
when there is strong cancelation between the S± , 5*2 , S3 
and S4 contributions to the reduced matrix element. For 
example, the largest difference is for the F% — 2 to F2 = 3 
transition in 229 Ra + . The data in Table fVl shows that 
the final value of the reduced matrix element for this 
case is just about 40% of the Si contribution. On the 
contrary, if the amplitude is dominated by the Si term 
the agreement between results of the two works is much 
better. This should be expected since the Si term is the 
most stable in the calculations. 



VI. CONCLUSION 

We present simultaneous calculation of the spin- 
independent and spin-dependent PNC amplitudes of the 
s-d transitions in Ba + , Yb + and Ra + . The results are 
to be used for accurate interpretation of future measure- 
ments in terms of the parameter of the spin-dependent 
PNC interaction x. Both, sign and value of x can be de- 
termined. Theoretical uncertainty is at the level of 3 to 
6% for Ba+ and Ra+ and 30 to 50% for Yb+. Note that 
the uncertainty for the spin-independent PNC amplitude 
can be further reduced by including structure radiation 
and ladder diagrams for more accurate treatment of cor- 
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relations and by including other small corrections (Breit, 
QED, etc.). The uncertainty for the relative contribution 
of the nuclear-spin-dependent part of the PNC amplitude 
is already small being on the level of 1% in some cases. 
The ratio of the SD to SI PNC amplitude is to be mea- 
sured to extract the calue of h. The results of the PNC 
calculations for Yb + are presented for the first time. 
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